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In ordinary statistical mechanics the Bohzmann-Shannon entropy is related to the Maxwell- 
Bolzmann distribution pi by means of a twofold link. The first link is differential and is offered by 
the Jaynes Maximum Entropy Principle. Indeed, the Maxwell-Boltzmann distribution is obtained by 
maximizing the Boltzmann-Shannon entropy under proper constraints. The second link is algebraic 
and imposes that both the entropy and the distribution must be expressed in terms of the same 
function in direct and inverse form. Indeed, the Maxwell-Boltzmann distribution pi is expressed 
in terms of the exponential function, while the Boltzmann-Shannon entropy is defined as the mean 
value of — In(pi). 

In generalized statistical mechanics the second link is customarily relaxed. Of course, the gener- 
alized exponential function defining the probability distribution function after inversion, produces 
a generalized logarithm A(pi). But, in general, the mean value of — A(pi) is not the entropy of the 
system. Here we reconsider the question first posed in [Phys. Rev. E 66, 056125 (2002) and 72, 
036108 (2005)], if and how is it possible to select generalized statistical theories in which the above 
mentioned twofold link between entropy and the distribution function continues to hold, such as in 
the case of ordinary statistical mechanics. 

Within this scenario, apart from the standard logarithmic-exponential functions that define or- 
dinary statistical mechanics, there emerge other new couples of direct-inverse functions, i.e. gener- 
alized logarithms A(a;) and generalized exponentials A~^{x), defining coherent and self-consistent 
generalized statistical theories. Interestingly, all these theories preserve the main features of ordinary 
statistical mechanics, and predict distribution functions presenting power-law tails. Furthermore, 
the obtained generalized entropies are both thermodynamically and Lesche stable. 

PACS numbers: 05.90.-|-m, 05.20.-y, 51.10.-|-y, 03.30.+P 



I. INTRODUCTION 

Apart from the exponential distribution exp(— a;), sev- 
eral distribution functions have been considered both 
from the theoretical point of view as well as to analyze 
experimental data. Ref. [1] lists more than 130 different 
distributions functions, more or less frequently used in 
statistical sciences. 

In the last decades particular attention has been de- 
voted to distribution functions exhibiting power-law tails, 
namely Ax'"" for x ^ oo, appearing in the phenomenol- 
ogy of many physical, natural and artificial systems. The 
power-law tailed distributions have been observed, for in- 
stance, in high energy physics (plasmas [2|, cosmic rays 
[1, 01 and particle production processes [aJMl; in con- 
densed matter physics (anomalous diffusion fluid 
motion llOl . turbulence 111 13), in natural sciences (seis- 
mology [lj|, meteorology [l4| and heliophysics [III), in 
biology (botany [T6j , genomics [l3l ) , in economics (stock 
prices [H, [l^ and personal annual income [131), in so- 
ciodemography (population of cities l2ll , frequencies of 
family names [23] , telephone calls [23] ) , in linguistics 
[13, [25'], in philology (classical mythology [13]), and also 
in fields involving human activity in general (intensity 
of wars [27I, terrorism attacks |2a, citation of scientific 



papers [29|, traffic in complex networks [30], etc.). 

The very simple analytic expression for the tails i.e. 



A'. 



in the above distributions, known as Pareto's or 
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Zipf's law, has been reg ularly adopted in the literature 
for data analysis [Slj, [sj, HI] . On the contrary, the ques- 
tion regarding the form of the analytic expression for the 
power-law tailed distribution that holds for < a; < +cxi 
remains, at the moment, an open problem. 

In the last few years there have been studies of the 
general properties of distribution functions of arbitrary 
forms [13, m, m, [13, HI]. On the other hand, the gen- 
eral properties of an arbitrary entropic form are also 
well known [13, 113, SH . Finally the connection bettween 
an arbitrary entropy and the relevant distribution given 
through the Jaynes maximum entropy principle |,42j , has 
been considered by several authors [13, 113| • 

A mechanism frequently used to explain the occurrence 
of statistical distributions is based on evolution equa- 
tions, essentially nonlinear Fokker-Plank but also Boltz- 
mann equations ^ HliiliE 111111113, 113, HH, 113, 
Issl . [13] ■ After fixing the form of the evolution equation 
the entropy of the system is automatically fixed and the 
probability distribution function is obtained as a station- 
ary state of the evolution equation. 

Clearly the correctness of an analytic expression for a 
given power-law tailed distribution, used to describe a 
statistical system, is strongly related to the validity of 
the generating mechanism. In the present paper we will 
consider power-law tailed distribution functions from a 
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general prospective and independently on the particular 
generating mechanisms. 

We recall that, according to Maximum Entropy Princi- 
ple, the Boltzmann-Shannon entropy S — — ^ p. ln(p. ), 
yields the exponential Maxwell-Boltzmann distribution. 
In generalized statistical mechanics, a special role is 
played by the trace-form entropies 



(1.1) 



where A{x) is an arbitrary, strictly increasing function 
defined for x > 0, which for x ^ 0^ becomes — oo, while 
p = {pi ,1 < i < N} is a discrete probability distribution. 
The function A{x) can be viewed as a generalization of 
the ordinary logarithmic function. We remark that the 
entropy is the ordinary mean value of —A{pi). It is clear 
that the form of A{x) imposes the form of the distribution 
function according to the maximum entropy principle. 

In proposing generalized statistical theories some stan- 
dard properties (positivity, continuity, concavity, symme- 
try, etc.) are customarily required for the entropy. Un- 
fortunately, these properties are not sufficiently strong to 
impose the form of the entropy; then, the metaphor of 
ordinary statistical mechanics remains the unique guid- 
ing principle for further generalization. In this sense the 
maximum entropy principle, the cornerstone of statisti- 
cal physics, is a valid and powerful tool to explore new 
roots in searching for generalized statistical theories. 

The main goal of the present paper is to show that the 
maximum entropy principle, by itself, suggests a gener- 
alization of ordinary statistical mechanics. The obtained 
generalized statistical theories preserve the main features 
of ordinary statistical mechanics and predict distribution 
functions, showing power-law tails, some of which are al- 
ready known in the literature. 

This paper is organized as follows: 

In Sect. 11 we consider the maximum entropy principle 
within the generalized statistical theories, based on trace 
form entropies. 

In Sect. Ill we obtain a class of entropies depending 
on three real parameters. 

In Sect. IV we consider the statistical distributions, re- 
lated to this class of three-parameter entropies, exhibit- 
ing power-law tails. 

In Sect. V and Sect. VI we discuss the thermodynamic 
stability and the Lesche stability of the generalized sta- 
tistical theories, based on the three-parameter entropies 
obtained in Sect. 111. 

In Sect. Vll we show that the class of the three- 
parameter entropies contains, as special cases, some al- 
ready known in the literature, entropies. 

In Sect. VIII we study a new, two-parameter subclass 
of entropies, belonging to the three-parameter class of 
entropies. 

In Sect. IX we make some concluding remarks. 

Finally in the Appendix we solve a differential- 
functional equation which permits us to obtain the three- 
parameter class of entropies. 



II. THE MAXIMUM ENTROPY PRINCIPLE 

Let us consider the generalized entropy defined in Eq. 
(|l.ip . We note that the particular probability distribu- 
tion given hy p = {Sia , 1 < i < N} where a is a fixed in- 
teger with 1 < a < TV, describes a state of the system for 
which we have the maximum information. For this state 
we set iS* = and this condition imposes for the general- 
ized logarithm (i) A(l) = and (ii) 0+ A(0+) = 0. Fur- 
thermore, in analogy with the ordinary logarithm we nor- 
malize the generalized logarithm through (iii) A'(l) = l. 

We introduce the constraints functional 



C — a^ 



1 



M-^ g{i)p. 



(2.1) 



where a„ and a — {a^, a^,..., ai} are the I + 1 La- 
grange multipliers, while g(i) = {g, (i), g^ (i), gi(i)} 
are the generator functions of the I moments M = 
{Mi,M2,...,M,}. 

The variational equation 



S 

Spi 



'S + C) =0 , 



(2.2) 



implies the maximization of the entropy S under the con- 
straints, imposing the conservation of the norm of p 



(2.3) 



and the a priori knowledge of the values of the I moments 
M = J2s{i)p^■ (2.4) 



Eq. (|2.2[) represents the maximum entropy principle and 
yields the relationship 



■S-PtA{pi) = -a - g(i) 

OPi 



(2.5) 



defining unambiguously the distribution function pi . 

We explain the meaning of the Lagrange multipliers &j 
by observing that a^ = aj(Mi, ... , M;) with j = 0, 1, ^ 
and S — S'mi,...,m, f38]. In the first step we multiply 
Eq. (j2.5p by Pi, then sum with respect to the index i, 
and finally we derive with respect to with /c = 1, /, 
obtaining 

^ ^ P. A{p,) + l^p,^^—p, A(k) 

, - , (2.6) 



dMk dpi 
dsLQ d 



dMk 
and then 



dMk 



dS 

dm. 



. P^dMkW.''' ^'''^' 



dao d 



I 



dMk dM, 



^a,M,-0 



(2.7) 
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In the second step we derive Eq. (|2.5p with respect to , 
subsequently multiply by pi and then sum with respect 
to the index i, obtaining 

Ed d , , , dan J-^ 9a, , 



dMk ^ dMk 



dMk 

Finally, by combining Eqs. ()2.7p and ()2.8p we have 



/ d d&i 
a, M, - — ^ M,- 



dMi 



' Yl g^, = - H 5jfc = -afc . (2.9) 



Then, for j = 1, Z it holds 

dS 



dM, 



(2.10) 



regardless on the form of the generalized entropy. 

III. THREE-PARAMETER ENTROPIES 

In the following, we are interested to the particular 
class of probability distribution functions which are given 
in terms of the generalized exponential function 



namely of the form 



a ■ g(i) - a, + 7? 
A 



(3.1) 



(3.2) 



a, A and rj being three scaling parameters. We note that 
the parameters A and r? produce a scaling to the Lagrange 
multipliers. 

Ordinary statistical mechanics corresponds to the 
choice P^{x) = ln(a;), £{x) = exp(a;) and {a = 1, A = 
1, 77 = 1}, or alternatively {a — 1/e, A = 1, 77 = 0}. We 
pose now the question regarding the possible existence of 
further couples of functions A(a;), £{x) and set of param- 
eters {a. A, 77} defining statistical theories, through Eqs. 
p.ip . (|2.5p and p.2p . different from ordinary statistical 
mechanics. 

Regarding this proposition we note that Eq. (|3.2p . after 
inversion, can be written in the form 

A A{pi/a) + r; = -a • + a^ , (3.3) 

which after comparison with Eq. (|2.5p yields 

d 



dp,, 



PiA{pi) = A A{pi/a) + 77 



(3.4) 



The latter equation is a first order differential- 
functional equation whose solutions determine unam- 
biguously the form of the generalized logarithm after 



taking into account the conditions (i) A(l) = 0, (ii) 
0+ A(0+) = 0+ and (iii) A'(l) = l. 

After obtaining the generalized logarithm and by inver- 
sion the generalized exponential, through Eqs. (|l.ip and 
p.2p . the entropy and the probability distribution func- 
tion are unambiguously fixed. 

Regarding the meaning of Eq. (|3.4p we note that the 
equation guarantees a twofold link between entropy and 
the probability distribution function: firstly we have a 
differential link offered by the maximum entropy prin- 
ciple through Eq. (|2.5P : secondly we have an algebraic 
link offered by Eq. (|l.ip and Eq. (|3.2p where the same 
function in direct and inverse form, i.e. A(a;) and £{x), 
defines the entropy and the probability distribution func- 
tion according to the standard rules of ordinary statisti- 
cal mechanics. Eq. p.4p has been proposed in ref. [s^ . 
while the special case corresponding to 77 = has been 
considered in ref. [5^ . The solutions of Eq. (|3.4I) for the 
special case 77 = 0, have been obtained in Ref. 36, 53]. 

Preliminarily we note that Eq. p.4p admits the clas- 
sical solution corresponding to A(a;) = ln(a;) and {a ~ 
1, A = 1, 77 = 1}, or alternatively {a — 1/e, X — 1, 77 = 
0}. The general solution of the Eq. p.4p . derived in the 
Appendix, is given in terms of a three-parameter function 
A{pi) — lnKr?(Pi) having the form 



ln^r^{x) 



(k + t)<J'' + {Hi - t)<^"« 



(3.5) 



The three parameters {k,t, <;} are related to the param- 
eters {a, A, 77} of Eq. p.4p through 



1 + r - K 
1 + r + K 



(1 



(1 + T + K 

77-(A-l) 



(k + r)?*^ + {k, — t)^ 



(3.6) 
(3.7) 
(3.8) 



The classical solution A{pi) — In(pi), a — 1/e, A = 1, 
7; = is contained as a limiting case (t = 0, k — > 0) in 
the above three parameter class of solutions. 

Let us consider the dual function of f{x) defined as 
follows: 



f(x) = -f{l/x) 



(3.9) 



It is easy to verify that the dual of a three-parameter 
logarithm again results to be a new three-parameter log- 
arithm, i.e. 



h\^r<;{x) = ln„^/^/(a;) , 



(3.10) 



where t' — —t and = l/<;. In general a three- 
parameter logarithm is not self-dual, but the class of the 
three-parameter logarithm is self-dual. 
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In the present context a special role is played by 
the scaling operation. Let us consider the class of the 
functions T — {f{x), x > 0}, satisfying the condition 
/(I) = and /'(I) — 1. Starting from the function 
f{x) G we define the corresponding scaled function 
<jf{x) G as follows: 



.f{x) 



f{<yx) - f{a) 



(3.11) 



where a > is the scaling parameter. 

It is remarkable that two consecutive scaling opera- 
tions are equivalent to a unique scaling operation, i.e. 



(3.12) 



where = (72 ci. 

The scaling operation has the following group proper- 
ties: 



i/(x) = f{x) , 



02 L o-i 



(3.13) 
(3.14) 
(3.15) 



It is straightforward to show that the three-parameter 
logarithm ln„T-t;(a;) has the interesting property 



■\nf^r<;{x) = \n^r^'{x) , 



(3.16) 



where ^' = cr<j. The scaled function cr^T^KTq{x) of a 
three-parameter logarithm again results to be a three- 
parameter logarithm with the same first two parameters, 
while the third parameter, given by = cr is the scal- 
ing parameter. In general a particular three-parameter 
logarithm is not self-scaling, but the class of the three- 
parameter logarithms h\KT<;{x) is self-scaling. 

The explicit form of the three-parameter class of en- 
tropies is given by 



S ^ - ^Pi hlKrc-CPi 



(3.17) 



IV. POWER-LAW TAILED DISTRIBUTIONS 

It is immediately possible to verify that In^r^Ca;) € 
C°° {R^) presents the following asymptotic power-law be- 
haviour: 



lnKT?(a;) 



A 



ln«,,(a;) ^ , 

X— i- + 00 



(4.1) 
(4.2) 



where A = [(|k| -I- t) c;!''! -I- (|k| - r) \ For r > 

|k| — 1 we have lim^_,Q+ x In^^Ttiix) = so that it results 



lnKr<;(a;) dx 



1 + (t^ - n") 



1 



(1 + r)2 - k2 



(k + t) i,'^ + [k. - t) 



(4.3) 



From Eqs. JSJ]) and ((37TO| it follows 

lnKr<;(l/a;) = -ln^^'^/(a;) , 



(4.4) 



whith r' = — T and q' = l/<j. 

After noting that when — |k| < r < |k| it holds 
d\nK,T^{x)ldx > 0, we can define the three-parameter ex- 
ponential exp^^^(a::) as the inverse function of \t1k,t^{x). 
The properties of exp^^^(a;) readily follow from those of 
\t1i^t^[x). For instance, the asymptotic power-law behav- 
ior of exp (x) is given by 



exp«^^(a;) - (±Ax) 



l/(r±|K|) 



X — y±oo 

while Eq. (|4.4|) implies 

exp^^^(a;) exp^^,^,(-x) = 1 



(4.5) 



(4.6) 



The probability distribution function can been ob- 
tained after maximization of the three-parameter entropy 
(I3.17p . according to maximum entropy principle 



a ■ g(0 ~&o+-n 
A 



(4.7) 



We remark that Eqs. dS^]), (|XT7| and (gTD define 
a very wide class of statistical theories, candidates to 
generalize the ordinary statistical mechanics after passing 
a series of validity tests, some of which will be considered 
in the next sections. 



V. THERMODYNAMIC STABILITY 

We denote with q = {qi\ the optimal distribution 
which according to maximum entropy principle, is de- 
fined by the algebraic equation (|2.5p . i.e. 



d 

qi lnKr<;('7j) = -a • g(z) + . 



(5.1) 



After introducing the entropy density for a general 
probability distribution function pi according to 

o'(Pi) = -Pi lnKr?(Pi) , (5.2) 
the constrained entropy $ = S* + C can be written as 



while its maximum value is given by 
da{qi 



dqi 



ao-a-M, (5.3) 



ao - a-M . (5.4) 



In order to calculate the difference of the above two func- 
tionals 



$(p)-$(<z)-E 



da{q,) 

<^\Pi)-(^\qi) ^ (Pi-ft) 

oqt 



(5.5) 
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when Pi « qi, we perform the Taylor expansion 



and obtain 



^(7iq^) 
dqi 



(pt - qi) + 



1 



{pi -qi? 



- $(g) = 



1 9Vfe) 



(pi - qi? 



(5.6) 



(5.7) 



At this point we recall that $(gi) represents the maxi- 
mum value of $(j3i) so that it holds 

$(p) - $(g) < . (5.8) 

The latter inequality and Eq. (|5.7p imply 

9V(x) 



< , 



(5.9) 



for any value of x. After taking into account the defini- 
tion of the density entropy, we get the following property 
for the generalized logarithm 



dx 



2 X \ri^r^{x) > 



(5.10) 



The latter relationship can be verified easily by starting 
with Eq. (|3.4p and after taking into account that for 
— |k| < r < |k| it results d\nnr';{x)/dx > 0. Therefore 
we can conclude that the system described by the entropy 
()3.17p is thcrmodynamically stable for — |k| < r < 



where we have used the definitions [x]+ = max(x, 0). The 
definition of A(p, s) implies that for any s the following 
relationship holds 



N 



A{p, S) - A{q, s)\<J2\P^-1^\ = \\P - 



(6.4) 



i=l 



while for s > —A lnK,-<;(l/A^) it results 

A(p, s)-A(g, s)| <iVaexp„^^(-^) . (6.5) 

After taking into account the definition of A{p, s) and 
Eq. p.4p , it is easy to verify that the entropy (|3.17p can 
be written in form 



-|-oo 

Sip)= J [1-A{p,s)]ds -1-7? 



(6.6) 



From the latter relationship it follows that the absolute 
difference of the entropies of the two different distribu- 
tions p — {Pi} and q — {q.} verifies the inequality 



Sip) - S{q) 



[A{p, s) - A{q, s)] ds 



< J \a{p, s)-A{q, s) 



ds 



VI. LESCHE STABILITY 

The Lesche stability condition [H, HO, |4l|, imposes 
that any physically meaningful entropy, depending on a 
probability distribution function, should exhibit a small 
relative error 



S{p) - S{q) 



(6.1) 



with regard to small changes in the probability distribu- 
tions 



D=\\p~q\ 



(6.2) 



Mathematically this implies that for any e > there ex- 
ists 5 > {) such that R < e holds for all the distribution 
functions satisfying D < S. It is known that the Lesche 
stability condition holds for the Boltzmann-Shannon en- 
tropy. 

In the following, by adopting the procedure of ref. [36l | , 
we will show that the Lesche stability condition holds for 
the entire class of three- parameter entropies (I3.17p . 

By using the procedure of ref. [s^ , we introduce the 
auxiliary function 



N 

A{p, s) = Y, [p. 

i=l 



a exp^ 



(-3 



(6.3) 



A{p, s) - A{q, s) ds . (6.7) 



Choosing -A lnKr?(l/A^) < ^ < +00, by using Eq fOl) 
in the first integral and Eq. (|6.5p in the second integral 
of Eq. (|6.7p . we obtain 



S{p)-S{q) < {£+!) lb -gill 

+ 00 



+ Na J exp^^^ ds . (6.. 



In particular the latter relationship holds for £ — £ being 

'l|p-9lli~ 



= — A In, 



aN 



(6.9) 



the values o{£ that minimizes the r.h.s. of the expression, 
as long as £ > — A lnKT<;(l/A^), which is true when 



S = lb - ^lli < a 



(6.10) 



i.e. for sufficiently close distributions, according to the 
metric. After some simple algebra one obtains the rela- 
tive difference of entropies 



S{p) - S{q) 



<ei6,N) 



(6.11) 



6 



and 
e{5, N) 



ln,,,(l/iV) 



In, 



■i6/N)-l + j 



(6.12) 



where S'max = ^ 1i1kt?(1/-^)- After taking into account 
that hm2,^o+ ^ liiKT?(a;) = 0, we find 



lim e(S,N) = 0+ 



(6.13) 



such that if the two distributions are sufficiently close, 
(5^0+, the corresponding difference e of entropies (|3.17|) 
can be made small at will and then we can conclude that 
the entropy (|3.17p is Lesche stable. 

In the thermodynamic limit N oo we have 



lim e{6,N) = e{S) = 6^- 

N—foc 



K +T 



(6.14) 



and results that lim5^o+ ^i^) = 0^. Thus the entropy 
(I3.17P is Lesche stable also in the thermodynamic limit. 



The properties of the functions ln„(x) and exp^(x) 
have been study extensively in the literature due to their 
relevance in relativistic statistical mechanics. The de- 
formation mechanism introduced by the parameter k is 
induced directly by the Lorentz transformations of Ein- 
stein's special relativity. Indeed 1/k is the reciprocal of 
light speed while 1/k^ is the particle rest energy, in di- 
mensionless forms. On the other hand the self-duality 
relation given by Eq. (|7.4|) leads to the relativistic dis- 
persion relation, obtaining in this way a direct link with 
the microscopic dynamics of the system. An important 
property of exp^(x) is given by 



exp„(9^ ©9b) = exp^(g^)exp^(gg) 



(7.6) 



1a ® 1b being the additivity law of the dimensionless 
relativistic momenta, defined through 



1a ffi 1b 



1 + 



(7.7) 



A further property of exp^(a;) is given by 



VII. SPECIAL CASES 

In the present section we consider a few special cases 
of the three parameter generalized logarithm p.Sp , some 
of which are already known in the literature. The ordi- 
nary logarithm, which is self-dual and self-scaling, can 
be obtained as the limiting case when k — > and r = 0, 
independently on the value of the parameter <j, such that 
we can write ln(x) — lnoo>;(a;). 

i) Self-dual logarithm: the choice corresponding to 
T = 0, = 1 and — 1 < K < 1, yields the generalized 
logarithm lnK(a;) = lnooi(a;), which together with its in- 
verse function assume the form 



Iuk^x) 



2k 



exp^(a:) = + k^x^ 



i/« 



(7.1) 
(7.2) 



exp«:(g) = exp^((7) 



where the K-derivative 



dq 



(7.8) 



(7.9) 



emerge within the Lorentz- invariant differential calculus. 

In order to better explain how the functions exp^(a;) 
and \ni^{x) can define a statistical theory we consider 
briefly the paradigm of statistical mechanics. We indicate 
with Wi the microscopic energy, with (3 the reciprocal of 
the temperature and with the chemical potential, so 
that, the Lagrange multipliers are given by ao = /3/x and 
ai = [3. The system entropy is given by 



S = - ^Vi InK(pj) , 



(7.10) 



The self-duality of the above functions is expressed 
through the relationships 



exp^(x)exp^(-a;) ^ 1 



(7.3) 
(7.4) 



It is remarkable that the self-duality property, 
if imposed to the three-parameter logarithm, i.e. 
liiKTi;(l/a;) = —\ni^T(:{^/x), automatically fixes the pa- 
rameters T = and <r = 1, so that h\K,{x) remains the 
only one self-dual generalized logarithm. 

The following integral representation for \ni^{x) holds 



ln„(x) 



1 



dt 



2 ./i/. 



(7.5) 



which in the limit of k — s- reduces to a well known 
property of the ordinary logarithm. 



while the maximum entropy principle yields the proba- 
bility distribution function 



a exp^ 



(7.11) 



fi being the normalization constant, while the expressions 
of the parameters A and a are given by Eqs (|3.6p and (|3.7p 
respectively, after posing r = 0. 

Statistical mechanics based on Eqs. (17.101) and (|7.1ip 
has been introduced in [sil, , while in |53l.l53| its rela- 
tivistic origin has been shown. 

In the last few years various authors have considered 
the foundations of this statistical theory, e.g., H-theorem 
and molecular chaos hypothesis 5^ 53, thermodynamic 
stability [Ullil, Lesche stabilitjTlsalea, [6ll, [il, Legen- 
dre structure of ensued thermodynamics [63[ jGibbs the- 
orem [gj], geometrical aspects of the theory [65| etc. On 
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the other hand specific apphcations to physical systems 
have been considered, e.g., cosmic rays [53], relativistic 
[66l | and classical [63| plasmas in the presence of external 
electromagnetic fields, relaxation in relativistic plasmas 
under wave-particle interactions (gsI . [69j , kinetics of in- 
teracting atoms and photons [70|. particle systems in ex- 
ternal conservative force fields [71l| , astrophysical systems 
[zl, [ill , quark-gluon plasma formation [74] , relativistic 



quantum hadrodynamics [75 
concern nonlinear diffusion 
the edge of chaos (t^ 



75 



etc. Other applications 
761 ] ■ dynamical systems at 
79], fractal systems [13], ran- 
dom matrix theory [81|, l82], error theory [s^, game the- 
ory [11], information theory [sl], etc. Also applications 
to economic systems have been considered, e.g., to study 
personal income distribution [86', 'st', ^8§\ and to model 
deterministic heterogeneity in tastes and product differ- 
entiation [H, [qO] etc. 

ii) Self-scaling logarithm: The case corresponding to 
T — —K — a/2, is undoubtedly the more discussed in 
the literature and produces the following generalized self- 
scaling logarithm and exponential: 



ln„(x) 



1-x- 



The self-scaling property is given by 
a-ln„(.T) = lnQ(a;) . 



(7.12) 
(7.13) 

(7.14) 



Both functions (|7.12p and (|7.13p have been considered 
first by Euler in 1779 and explicitly has remarked that the 
function {x'^ — 1)/lu reduces to the natural logarithm in 
the w — > limit [9l|. Successively, the function (|7.13[) 
was adopted in 1908 by Gosset to construct the Stu- 
dent distribution [93], which in 1968 was introduced into 
plasma physics by Vasyliunas Q . In mathematical statis- 
tics the function (17.131) has been used to construct the 
Burr-distribution [93]. On the other hand the function 
(|71^ was adopted in 1967 by Harvda e Charvat [9^] to 
propose a generalized entropy in information theory. Fi- 
nally both the functions (|77T2l) and fTTS)) were adopted 
by Tsallis in 1988, using the new parameter q = a + 1 to 
define the q- logarithm hiq{x) and q-exponential exp^(a;). 



\llq{x) 



c^-l - 1 
1-g 



(7.15) 



exp^(x) = (l + (l-g)x)i/(i-') , (7.16) 
in order to develop nonextensive statistical mechanics 

[IE pi]. 

iii) The case corresponding to the choice k = (a — 
l/a)/2, r = (a + l/a)/2 — 1 and <; = 1 defines a gener- 
alized logarithm 



lna(a;) 



^l/a-l 



(7.17) 



a — 1/a 

leading, by means of Eq. (I3.17|) . to an entropy intro- 
duced in the hterature in 1997 by Abe [97j] presenting 
the symmetry a <-> 1/a. 



iv) Only in very few cases, the generalized logarithm 
given by Eq. p.5p , can be inverted to obtain analytically 
the corresponding generalized exponential. For instance 
the choice k — 37/2, r — 7/2, = 1, yields an invertible 
generalized logarithm 



In (x) 



r27 



37 



corresponding to the generalized exponential 



exp (x) 



1 + y 



1/3 



1-2/ 



1/3- 



.1 1/7 



(7.18) 



(7.19) 



with y = v/1 - 473a;3, [sg. 

v) The three-parameter class of generalized logarithm 
(13. 5p contains a two-parameter subclass, already known 
in the literature, defined by posing ^ — 1 



hlf^rlix) = X'^ IUk^x) 



(7.20) 



The latter subclass of two-parameter-logarithms is not 
self-scaling. From Eq. (|3.16p with <j = 1/^, one obtains 
the following relationships: 

ln„^<;(x) = ^Irif^riix) , (7-21) 
ln„^i(a;) = i/^lii,^r<;{x) . (7.22) 

After introducing the new parameters a and (3 through 
K = (a — /3)/2 and t = [a + /3)/2, the two-parameter 
logarithm hif^Ti{x) = hia(i{x) can be written in the form 



In Jx) 



a- (3 



(7.23) 



which is more familiar in the literature and has been 
considered firstly by Euler [91]. In general this two- 
parameter logarithm cannot be inverted analytically, 
thus it is impossible to define explicitly the correspond- 
ing generalized exponential. The function (|7.23[) has 
been used to construct the two-parameter entropies in- 
troduced in the information theory in 1975 by Mittal 
[9^, and by Sharma and Taneja [9^. The generalized 
logarithm (|7.23p was reconsidered i n st atistical mechan- 
ics in 1998 by Borges and Roditi [lOOl ] and recently in 
refs. |35j, [3^. Clearly the k- logarithm does not emerge 
easily from (|7.23p and probably this is the reason why it 
has not been discussed in the literature before 2001. 

vi) The three-parameter class of logarithm (|3.5|1 con- 
tains a second, new two-parameter subclass of loga- 
rithms, namely the scaled K-logarithms, obtained by pos- 
ing T — 0. The corresponding class of generalized en- 
tropies and distribution functions will be discussed in the 
next section. 

The one-parameter and two-parameter generalized log- 
arithms, already known in the literature and discussed in 
the present section, emerging as special cases of the three- 
parameter logarithm (13. 5p . are reported in the table. 



8 



TABLE L Generalized logarithms already known in the literature, obtained here as special cases of the three-parameter 
logarithm. In the third column are reported the values of the parametes n, r, and c to insert in the general expression of the 
three-parameter logarithm, to obtain as special cases the zero-, one-, and two-parameter logarithms. 



generalized logarithm expression 



parameter values 



introduced in ref. 



three-parameter logarithm 
(most general case) 



InKr^(x) 



(k + t)<;'= + {k-t)5-'« 



1. zero-parameter logarithm 



ln(a;) 



r = 0, K-»0 



Napier 



2. one-parameter logarithm 



ln4x) = si 



r = 0, ? = 1 



[51J 



3. one-parameter logarithm 



K = T = {l-q)/2 



[26] 



4. one-parameter logarithm 



Inafa;) 



X —X 



a — l/a 



? = K=(l-l/a)/2 
r = (l + l/a)/2-l 



[97] 



5. two-parameter logarithm 



T = {a + l3)/2 



[98, 99] 



6. two-parameter logarithm In^zix) = ^—^ — ^^^^^^„^j 



r = 



VIII. SCALED «:-STATISTICAL MECHANICS 

We consider here the two-parameter generalized loga- 
rithms hiK^(x) = In K 0^(2;), obtained by posing t = in 
Eq. (|3.5[) . It is easy to verify that lni^^{x) is the scaled 
function of In^ (x) . The following relationships hold 



lnK?(a;) = <; ln„(x) , 
lnK(.T) = i/^\n,,^{x) , 
ln,^{x) = hi„i(x) , 



(8.1) 
(8.2) 
(8.3) 



and results lno;;(a;) — ln(a;). The scaled k- logarithm as- 
sumes the form 



lnK?(a;) 



and can also be written as 



ln„<;(x) = - ln„(<^a;) 
a 



(8.4) 



(8.5) 



where a = -y^l -I- K-^[lnK(<;)]^ and b = InK(^). An impor- 
tant property of the scaled K-logarithm is that its inverse 
function, the scaled K-exponential, exists and assumes a 
very simple form for any value of k and <;: 



exPK^(a;) = - exp^(aa; + b) 
The following Taylor expansions hold: 



1 - - — 



ln^^(l -I- a;) x 



exp^^(x) ^ 1 + a; + ( 1 - - ^ 



(8.6) 



■7) 



(8.8) 



while the power-law asymptotic behaviours of ln„<;(.T) 
and exp^^(a;) are given by 



InK^(x) 



X 



lnK5(a;) 



X — > + oo 



2a|K| 

2,|k| 

2aH 



(8.9) 
(8.10) 
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exp«,(a;) 



\2aK 



±i/kl 



11) 



The new class of two-parameter entropies is given by 
S = -Y,P^ ln«,(p.) . (8.12) 

i 

Consequently the maximum entropy principle produces 
a distribution function having the form 



Pi= a exp„ 



a ■ g(^) - a„ + r/ 
A 



with 



1 - K 



1/2k 



»7 = 



1 - r 



(8.13) 

(8.14) 
(8.15) 
(8.16) 



The K-logarithm results to be the semi-sum of two 
scaled ^-logarithms with scaling parameters and 1/^ 
respectively, i.e. 



\n^{x) = i [ln„.;(a;) + ln«, 1/5(0;)] 



.17) 



A further connection between the K-logarithm and the 
scaled K-logarithm is given by 



ln„<;(a;) 



1 



\n,{x) + -\n^,,{x) , (8.18) 



with 1;' = 1/a. The following two expressions of \ni^^i{x) 
hold 

ln«,,,(a;) = Aln«(x/a) - 1 , (8.19) 
\Ti^,'{x)^ ^[x\ri^{x)]-l . (8.20) 

It is remarkable that \ni^^{x) in the limits of <; — > and 
— > cx) reduces to the self-scaling g— logarithm given by 

Ing(x) = (1 — x'^^'') /{q — 1). After posing q = 1 + |k| we 

obtain 



ln,(x) ^\iii^o{x) , 
ln2_,(x) = In K 00 (a;) 



5.21) 
5.22) 



On the other hand Eqs. ([8T7| and (|8?T8l) in the ^ ^ 
limit, reduce to 



\n^{x) = i [ lnq{x) +ln2-g(a;)" 



(8.23) 



ln,(.x) = 1 + ln,(x) - — ln,,,(x) , (8.24) 



with (7 — 1 + |k| and = 1/a. 



IX. CONCLUSIONS 

As conclusions, we recall briefly the main results ob- 
tained here. The requirement that the generalized trace- 
form entropy 



(9.1) 



after maximization, according to maximum entropy prin- 
ciple, yields a probability distribution function given by 



Pi = aA' 



1 f a ■ - a„ + 7] 
A 



(9.2) 



is sufficient to fix the form of A(a;) and therefore the 
form of the entropy and of the probability distribution 
function. 

The obtained class of generalized logarithms, A (a;) — 
^'0-KT^{x), depends on three free parameters according to 



In nT.;{x) 



{k + t)<;'' + (k - t)^"'^ 



(9.3) 



The relevant three-parameter entropy is both thermo- 
dynamically stable and Lesche stable and contains as spe- 
cial cases all the one-parameter and two-parameter trace 
form entropies appeared in the literature. 

The ensuing three-parameter probability distribution 
functions represent the minimal deformation of the 
Maxwell-Boltzmann exponential distribution compatible 
with the maximum entropy principle and exhibit power- 
law tails. 



X. APPENDIX 



Here we solve the equation 



d 



— xA(x) = A A(x/a) + rj 
ox 



(10.1) 



with the conditions A(l) = 0, A'(l) = 1, 0+ A(0+) = 
and A(0+) < 0. 

After expressing the function A{x) in terms of the aux- 
iliary function L{x) according to 



A{x) 



Lj^x) - L(^) 



and after setting 



(10.2) 



(10.3) 



we obtain the following equation for the function L{x) 
■a;L(a;) = Ai(^) . (10.4) 



d_ 

dx 
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The conditions associated to the function L{x) easily fol- 
low from those related to the function A(a;), i.e. L{1) = 0, 
L'(l) 1, 0+L(0+) = and L(0+) < 0. 

The general solution of Eq. (|10.4p has been obtained 
in Ref. [s^. In the following we recall briefly the main 
steps of the solution procedure. 

First we perform the change of variable L{x) = 
{l/x)h{Xa Inx) and x = exp(i/(Aa)) and reduce Eq. 
(|10.4p to a first order homogeneous differential-difference 
equation, belonging to the class of delay equations 



dh{t) 
dt 



(10.5) 



where to = Xa Ina. The general solution of this equation 
is given by 



i=l j=0 



(10.6) 



where n is the number of real solutions s . (of multiplicity 
mj of the characteristic equation s. — e^*o ''i = 0. In 
general the integration constants a. . depend on the 
parameter Si. The latter algebraic equation docs not 
admit real solutions for tg < — 1/e. For t„ > — 1/e 
the number and multiplicity of the real solutions of 
the characteristic equation depend on the value of to 
according to 

(a) for — 1/e < < 0, n = 2 and m. ~ 1, 

(b) for = —1/e, ri = 1 and m = 2, 

(c) for tg > 0, n = 1 and m = 1. 



The solution of p0.4p for the case (a) can be written 
in the form 

L(x) = Aj(k^, kJo;'"! + A,(Kj, Kjx'^a , (10.7) 

where k. = Xa s. ^ 1 and A. (k^ , ) the integration con- 
stants. The boundary conditions L{1) — and L'(l) = 1 



imply = —A^ and A^ + A^ = (k^ — k^) A-^ =1 
respectively. After introducing the new parameters k — 
(Kj ~k^)/2 and r = (k^ -^n^jl^ we can write the solution 
as follows 



L{x) 



2k 



(10.8) 



The condition L{Q~^) < implies r < From the 

characteristic equations we obtain the system 1 ± r -|- 
K = Xa^'^^'^, which can be solved to determine the two 
constants a and A. The constant 77 follows from Eqs. 
p0:3l) and ((Toll). We obtain 



a 



1 



(1 



(1+T 

(A 



(10.9) 



(10.10) 



1 



(k + r)<^« (k - r)c;-« ^ ^ 
with r > |k| — 1. This latter condition guarantees that 
0+L(0+) = 

Case (b) corresponds to the limit of case (a) 
producing L{x) — x'^ lux. Case (c) gives the trivial 
solution L{x) = ax^ with 6 = A a s — 1 which can not be 
used to define a generalized logarithm. 

The most general expression of L(x) is given by 
Eg. pas)) . Finally, after taking into account Eq. p^?2l) 
we can write the explicit form of the generalized loga- 
rithm A(a;) depending on the three parameters {k, r, <;}: 



A{x) 



(k + t)^'' + (k - t)^ 



(10.12) 



The conditions A(0+) < and 0+A(0+) = imply that 
T < |k| and r > — 1 respectively. 
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